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ABSTRACT 

In this paper, an M/M/l queueing system with server Start-up, N-Policy, breakdowns at two phases of service 
and Balking is analyzed. Customers arrive at the system according to a Poisson process. They receive the first 
phase service as a batch followed by second essential phase of individual service. Arriving customers may balk 
with a certain probability and may depart without getting service due to impatience. Lack of service occurs 
when the server is on vacation, busy during the service mode or due to the sudden breakdowns of the server. 
For this model the probability generating functions for the number of customers present in the system at various 
states of the server are derived and obtained the closed-form expressions for various performance measures of 
interest. Further a total expected cost model is formulated to determine the optimal threshold ofNat a minimum 
cost. Finally, numerical examples are given. 
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I. INTRODUCTION 

We consider the optimal policy of a removable and un-reliable server for an M/M/l two-phase 
queueing system, where the removable server operates an N policy. 

Two-phase queuing systems have been discussed in the past for their applications in various areas, such 
as computer, communication, manufacturing, and other stochastic systems. In many computer and 
communication service systems, the situation in which arriving packets receive batch mode service in the first 
phase followed by individual services in the second phase is common. As related literature we should mention 
some papers [1,7,12] arising from distributed system control where all customers receives batch mode service in 
the first phase followed by individual service in the second phase. 

The concept of the N policy, was first introduced by Yadin and Naor [13]. Past work regarding 
queueing systems under the N policy may be divided into two categories: (i) cases with server vacations, and (ii) 
cases with server breakdowns. 

Queueing systems with server vacations have attracted much attention from numerous researchers since 
the paper was presented by Levy and Yechiali[9]. Server vacations are useful for the system where the server 
wants to utilize his idle time for different purposes. An excellent survey of queueing systems with server 
vacations can be found in papers by Doshi and Takagi [3,4,5, 11]. Queueing models with server vacations 
accommodate the real-world situations more closely. Such model frequently occurs in areas of computer and 
communications, or manufacturing systems. For example, consider an assembly line where a worker may have 
some idle time between subsequent jobs. To utilize the time effectively, managers can assign secondary jobs to 
the worker. However, it is important that the worker must return to do his primary jobs when he completes the 
secondary jobs. 

Queueing systems with impatient customers have been studied by a number of authors [2, 10] .There is 
an extensive amount of literature based on this kind of model .The source of impatience has always been taken 
to be either a long wait already experienced upon arrival at a queue, or a long wait anticipated by a customer 
upon arrival. Haight [6] first considered an M/M/l queue with balking. 
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However, to the best of our knowledge, for two -phase queueing systems with N-Policy, server 
breakdowns, there is no literature which takes customers' impatience into consideration. This motivates us to 
study a two-phase queueing system with N-policy, server start-up, breakdowns and balking. Thus, in this present 
paper, we consider two-phase M/M/l queueing system with server Start-up, N-Policy, unreliable server and 
Balking where customers become impatient when the server is unavailable. 

The article is organized as follows. A full description of the model is given in Section. 2. The steady- 
state analysis of the system state probabilities is performed through the generating in Section. 3 while some, 
very useful for the analysis, results on the expected number of customers in different states and some special 
cases are given in Section. 4. In Section. 5 the characteristic features of the system are investigated. Optimal 
control policy is explained in section.6, while, in Section. 7, numerical results are obtained and used to compare 
system performance under various changes of the parameters through sensitivity analysis. 

The main objectives of the analysis carried out in this paper for the optimal control policy are: 

i. To establish the steady state equations and obtain the steady state probability distribution of the number of 
customers in the system in each state. 

ii. To derive expressions for the expected number of customers in the system when the server is in vacation, in 
startup, in batch service (working and broken conditions) and in individual service (working and broken 
conditions) respectively. 

iii. To formulate the total expected cost functions for the system, and determine the optimal value of the control 
parameter N. 

iv. To carry out sensitivity analysis on the optimal value of N and the minimum expected cost for various 
system parameters through numerical experiments. 

II. THE SYSTEM AND ASSUMPTIONS 

We consider the M/M/l queueing system with N-policy, two phases of service, server Breakdowns and 
balking with the following assumptions: 

1. Customers are assumed to arrive according to Poisson process with mean arrival rate X and join the batch 
queue. Customers will get the service in the order in which they arrive. 

2. The service is in two phases. The first phase of service is batch service to all customers waiting in the 
queue. On completion of batch service, the server immediately proceeds to the second phase to serve all 
customers in the batch individually. Batch service time is assumed to follow exponential distribution with 
mean 1//? which is independent of batch size. Individual service times are assumed to be exponentially 
distributed with mean On completion of individual service, the server returns to the batch queue to 
serve the customers who have arrived. If the customers are waiting, the server starts the batch service 
followed by individual service to each customer in the batch. If no customer is waiting the server takes a 
vacation. 

3. Whenever the system becomes empty, the server is turned off. As soon as the total number of arrivals in the 
queue reaches or exceeds the pre-determined threshold N, the server is turned on and is temporarily 
unavailable for the waiting customers. The server needs a startup time which follows an exponential 
distribution with mean 1/9. As soon as the server finishes startup, it starts serving the first phase of waiting 
customers. 

4. The customers who arrive during the batch service are also allowed to join the batch queue which is in 
service. 

5. The breakdowns are generated by an exogenous Poisson process with rates t>\ for the first phase of service 
and i for the second phase of service. When the server fails it is immediately repaired at a repair rate^ 2 in 
first phase and a 2 in second phase, where the repair times are exponentially distributed. After repair the 
server immediately resumes the concerned service. 

6. A customer may balk from the queue station with probability b 0 when the server is in vacation or may balk 
with a probability bi when the server is in service made due to impatience. 

III. STEADY-STATE ANALYSES 

In steady - state the following notations are used. 

Fb.L.o = The probability that there are i customers in the batch queue when the server is on vacation, where 
i = 0,1,2,3, ...,N-1 
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Pi.!.: = The probability that there are i customers in the batch queue when the server is doing pre-service 
(startup work), where i = N,N+l,N+2,. . . 

f3^ i [,o= The probability that there are i customers in the batch queue when the server is in batch service where i 



^,[,0= The probability that there are i customers in batch queue when the server is working but found to be 
broken down, where i = 1,2,3,... 

p4,i l j= The probability that there are i customers in the batch queue and j customers in individual queue when 
the server is in individual service, where i=0,l,2... and j= 1,2,3,... 

p^ . j -The probability that there are i customers in the batch queue and j customers in individual queue when 
the server is working but found to be broken down, where i = 0,1,2,. . . and j = 1,2,3, . . . 

The steady- state equations governing the system size probabilities are as follows: 
^oftAo = (1) 
- ; 1 < t < JY - 1. (2) 

(A6 L + e}p UN # = A6 0 ft, iJV _ li0 . (3) 



Ui L -h£ -f Oft.iio = ^ i} iP2,i-uo + M Pi.i.i + ^ 1 < t < JV -1- (5) 

Ui L + P + ^ L 3ft,i;o = ^ifr.i-uj + p P 4f l p l + frft ,i,o + 8Puu* 1 ^ N • (6) 

Ui L +■ a L +■ = -^frj.o + a 2p5,o,r; ^ l-(8) 

Ui> L + fl!+ /Op 4l E,j = m.ij+i +■ ^lP4,E-W + a 2?W E *J" ^ 1(9) 
Ui L + a 2 ^^j = fliP^J ^1- (10) 
Ui L -h a^ft.ij = fliP^uj + t_ 'i - 1- (H) 

To obtain the analytical closed expression ofpo lfll o, the technique of probability generating function can 

be successfully applied as detailed below. Define probability generating functions associated with marginal 
queue size distributions as follows: 

Co GO = SiLo 1 ^,^ 1 , f^GO = lf =N p Ll L lfl z L , 

Multiplication of equation (2) by z 1 and adding over i (1 < i < N-l) gives 



= 12 3 



01*! +■ £)p!,i,o = ^iPi,;-!^ t > AT- 



(4) 



(12) 
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Multiplication of equations (3) and (1.4) by z 1 and adding over i (i > N) gives 



(13) 



Multiplication of equations (1.5) and (1.6) by z 1 and adding over i (i > 1) gives 
(Ai L (l - z) + ^ +■ y C2 &>=^ GO +■ ^ ^GO +■ 0 ff t fc) - ^oftuui • 

Multiplication of equation (1.7) by z 1 and adding over i (i > 1) gives 



(14) 



(15) 



Multiplication of equations (1.8) and (1.9) by z L y J and adding over Corresponding values of i and j gives 



(16) 



Multiplication of equations (10) and (11) by z L y J and adding over Corresponding values of i and j gives 



Ui L (l - z) -h cc^G^C^y) = ■a 1 G 4 (z,y}. 

The total probability generating function G (z, y) is given by 
Giz> y) — G § fe) -h f7i Cz) -f -f f7 3 (z) -f G 4 (z, y) -f G 5 (z,y}. 

The normalizing condition is 

ffCU] = G 0 (1) + ff t (13 + Ga (1) -h C7 3 CO -f C 4 Cl,l] +■ GsCl.l) = 1 
From equations (12) to (19) 



ff 4 (U] = 



( 



■■2 -, (jfG± (LJ-.liflitL?)' 



let t L = _ 2 — Ai L (a L -f , then 



_____ 1+ P__J — ^ A< 



(17) 



(18) 



(19) 



C 0 CO = JVpM* , 


(20) 




(21) 




(22) 


c a C0 = (jj)fi*(D 


(23) 



(24) 



(25) 
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Substituting the value of Ri(l) from (22) to (25) gives G 2 (l), G 3 (1), G 4 (1,1) and G 5 (l,l). 

Probability that the server is neither in batch service nor in individual service is given by 
£*(1) + ClCD = 1 - f^(l + + a 

1 ^ P ^ W ft ^ «i 
This gives a, 0 n = (1 — p] . (26) 

Where p = (l -f -f ^ (l -f ^J^J is the utilizing factor of the system. 

From Equation (26) we have p < 1, which is the necessary and sufficient condition under which steady state 
solution exits. 

Under steady state conditions, let Pb.Fi.J^.J^.Fi, ond p 5 be the probabilities that the server is in vacation, startup, 

in batch service, in batch service with break down, in individual service and in individual service with break- 
down states respectively. Then, 

ft=ff 0 ClL (27) 

pi = ^ai (28) 

(29) 
(30) 

jj 4 = (31) 
Ps = ^(1,1). (32) 



4.1 Expected number of customers at different states of the server 

Using the probability generating functions expected number of customers in the system at different 
states are presented below. 

Let L 0 , Li, L 2 , L 3 L 4 and L 5 be the expected number of customers in the system when the server is in 
idle, startup, batch service, break down in batch service, individual service and break down in individual states 
respectively. 

Then 

= S^ 1 " tf^ l [ 1 o=^(l)= _ ^ i Fb,^o- (33) 

*i = ^l=jv t Pu^=^l(1)= f 2 FM*- (34) 

£ a -^l=l .uo=t' 2 (1)4 TVt J^ a 2- (35) 

£ a -^l=l ! == c '3(l)= ^ ■ (36) 
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2tn 



(37) 



=^%+^C 5 f 1,0.(38) 

The expected number of customers in the system is given by 
£(A0 = £ 0 + £i +£ 2 +£ a + £ 4 + £ 5 . (39) 

4.2 Special cases 

In this section we present some existing results in the literature which are special cases of our system. 

(i) When there is no startup, breakdowns in both the phases of service, balking and N=l,the equation (39) 
becomes 

(ii) £CftO + . (39a) 

This is the expected system length for the two-phase M/M/l queueing system without gating. The obtained 
expression coincides with that of Krishna and Lee (1990). 

(iii) When there is no batch service, the equation (39 a) becomes 
itiO = ^ (39b) 

This is the expected system length for the M/M/l queueing system. 

IV. CHARACTERISTIC FEATURES OF THE SYSTEM 

In this section, we obtain the expected system length when the server is in different states. Let 
E^ f E v E 2f E^ f E 4 and iT 5 denote the expected length of vacation period, startup period, batch service period , 

batch service breakdown period, individual service period, and breakdown period during individual service 
respectively. Then the expected length of a busy cycle is given by 

E c = E$ + F L + F 2 + F a + F 4 + F 5 . 

The long run fractions of time the server is in different states are as follows: 

r = n» ( 4 °) 

|=Pi, (41) 
jj=ft, (42) 

f^Pa, (44) 

Er 

t=K- (45) 

Expected length of vacation period is given by 
E ° = it W 
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Hence, 

E * = (JftaM ' (47) 

V. OPTIMAL CONTROL POLICY 

In this section, we determine the optimal value of N that minimizes the long run average cost of two- phase 
M/M/l, N-policy queue with server break downs and balking. To determine the optimal value of N, we consider 
the following cost structure: 

Let T (N) be the average cost per unit of time, then 

tcbo = c k m + c.f± + $ + c^) + c bL (D + c. 0 

(48) 
Where 

C k = Holding cost per unit time for each customer present in the system, 
C 5 = Cost per unit time for keeping the server on and in operation, 
Cp. = Startup cost per unit time, 
C a = Setup cost per cycle, 

C fil = Break down cost per unit time for the unavailable server in batch service mode, 

Cfo = Break down cost per unit time for the unavailable server in individual service mode, 
C b = Cost per unit time when a customer balks, 

C v = Reward per unit time as the server is doing secondary work in vacation. 

For the determination of the optimal operating N-policy, minimize T (N) in equation 48. 
An approximate value of the optimal threshold N* can be found by solving the equation 



dT 1 (N) 



dN 



(49) 

= o v y 



N = N * 



A computational algorithm translated in MATLAB is used to obtain the optimum values. 

VI. SENSITIVITY ANALYSIS 

In order to verify the efficiency of our analytical results, we perform numerical experiment by using 
MATLAB. The variations of different parameters (both monetary and non -monetary) on the optimal threshold 
N*, mean number of jobs in the system and minimum expected cost are shown. 

Parameters for which the model is relatively sensitive would require more attention of researchers, as 
compared to the parameters for which the model is relatively insensitive or less sensitive. 

We perform the sensitivity analysis by fixing 

Non -monetary parameters as 

^=0.5, Li=8,ai=0.2,a 2 =3. 0,^i=0. 2,^2=0.3, 9=6,p=12, b 0 =0.4 , bi=0.2 and monetary parameters as 
C r =15,C b i=50,C b2 =75,C b =15,C m =200,C h =5 and C s =1000; 
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7.1. Effect of variation in the non-monetary parameters 

Table 1: Effect of X on N*, expected system length and minimum expected cost 





0.9 


3 


1.7 


2.1 


2.5 


2.9 


N* 


12 


14 


16 


18 


19 


20 


L(N*) 


6 


7 


7 


8 


9 


9 


T(N*) 


51.58 


69.7 


81.97 


96.12 


107.39 


117.67 



It is observed from the above table that with increase in values of X, 

a) N* is increasing. 

b) Mean number of customers in the system is increasing. 

c) Minimum expected cost is increasing. 

Table 2: Effect of on \i andp N*, expected system length and minimum expected cost 





9,12.05 


10,12.1 
0 


11,12.1 
5 


12.12.2 
0 


13,12.2 
5 


14,12.3 
0 


N* 


9,9 


9,9 


9,9 


9,9 


9,9 


9,9 


L(N*) 


4,4 


4,4 


4,4 


4,4 


4,4 


4,4 


T(N*) 


31.71, 
31.92 


31.58, 
31.91 


31.48, 
31.91 


31.39, 
31.91 


31.31, 
31.91 


31.25, 
31.91 



It is observed from the above Table that, with increase in values of \i and p, 

a) N* is almost insensitive. 

b) Mean number of customers in the system is almost insensitive. 

c) Minimum expected cost is decreasing. 

Table 3: Effect of ai and a 2 on N*, expected system length and minimum expected cost 





0.205, 
1.1 


0.21, 
1.2 


0.215, 
1.3 


0.220, 
1.4 


0.225, 
1.5 


0.230, 
1.6 


N* 


9,9 


9,9 


9,9 


9,9 


9,9 


9,9 


L(N*) 


4,4 


4,4 


4,4 


4,4 


4,4 


4,4 


T(N*) 


31.92, 
31.87 


31.92, 
31.87 


31.92, 
31.86 


31.93, 
31.86 


31.93, 
31.86 


31.93, 
31.85 



It is observed from the above table that with increase in values of ai and a 2 , 

a) N* is almost insensitive. 

b) Mean number of customers in the system is insensitive. 

c) Minimum expected cost is slightly increasing withaiandslightly decreasing with a 2 . 



Table 4: Effect of 9 on N*, expected system length and minimum expected cost 



9 


7 


8 


9 


10 


11 


12 


N* 


9 


9 


9 


9 


9 


9 


L(N*) 


4 


4 


4 


4 


4 


4 


T(N*) 


31.72 


31.64 


31.58 


31.54 


31.51 


30.89 



It is observed from the above table that with the increase in values of 9, 

a) N* is decreasing. 

b) Mean number of customers in the system is decreasing. 

c) Minimum expected cost is decreasing. 
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Table 5: Effect of and ^ 2 on N*, expected system length and minimum expected cost 



^2 


0.3, 
0.305 


0.4, 
0.31 


0.5, 
0.315 


0.6, 
0.320 


0.7, 
0.325 


0.8, 
0.330 


N* 


9,9 


9,9 


9,9 


9,9 


9,9 


9.9 


L(N*) 


4,4 


4,4 


4,4 


4,4 


4,4 


4,4 


T(N*) 


31.87, 
31.92 


31.91, 
31.91 


34.00, 
31.91 


34.16, 
31.91 


34.39, 
31.91 


34.7, 
31.91 



It is observed from the above Table that, with increase in values of £i and § 2 ? 

a) N* is almost insensitive. 

b) Mean number of customers in the system is almost insensitive. 

c) Minimum expected cost is increasing with the increase of %\ and slightly decreasing with the increase of ^ 2 • 

Table 6: Effect of b 0 and bi on N*, expected system length and minimum expected cost 



b 0 , bi 


0.4015, 
0.25 


0.4030, 
0.30 


0.4045, 
0.35 


0.4060, 
0.40 


0.4075, 
0.45 


0.4090, 
0.50 


N* 


9,9 


9,9 


9,9 


9,9 


9,9 


9,9 


L(N*) 


4,4 


4,4 


4,4 


4,4 


4,4 


4,4 


T(N*) 


31.99, 
34.05 


34.07, 
34.17 


34.14, 
34.29 


34.22, 
34.41 


34.3, 
34.53 


34.37, 
34.64 



It is observed from the above Table that with the increase in values of b 0 and bi, 

a) N* is slightly increasing. 

b) Mean number of customers in the system is slightly increasing. 

c) Minimum expected cost is increasing with both b 0 and bi 

7.2. Effect of variation in the monetary parameters 

Table 7: Effect of C r on N*, expected system length and minimum expected cost 



Cr 


17 


19 


21 


23 


25 


27 


N* 


9 


9 


9 


9 


9 


9 


L(N*) 


4 


4 


4 


4 


4 


4 


T(N*) 


31.98 


30.04 


28.10 


26.16 


24.23 


22.29 



It is observed from the above table that with the increase in values of C r 

a) N* is slightly increasing. 

b) Mean number of customers in the system is slightly increasing. 

c) Minimum expected cost is decreasing. 

Table 8: Effect of Cbiand Cb2 on N*, expected system length and minimum expected cost 



Cbl, Q)2 


52,80 


54,85 


56,90 


58,95 


60.100 


62,105 


N* 


9,9 


9,9 


9,9 


9,9 


9,9 


9,9 


L(N*) 


4,4 


4,4 


4,4 


4,4 


4,4 


4,4 


T(N*) 


31.92, 
31.9 


31.92, 
31.94 


31.92, 
31.94 


31.92, 
31.95 


31.92, 
31.96 


31.92, 
31.97 



It is observed from the above table that with the increase in values of Cbiand Cb2, 

a) N* is almost insensitive. 

b) Mean number of customers in the system is almost insensitive. 

c) Minimum expected cost is insensitive with C b i and slightly increasing with C b2 . 
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Table 9: Effect of Cb on N*, expected system length and minimum expected cost 



c b 


20 


25 


30 


35 


40 


45 


N* 


9 


9 


9 


9 


9 


9 


L(N*) 


4 


4 


4 


4 


4 


4 


T(N*) 


35.43 


36.94 


38.463 


39.97 


41.49 


43.01 



It is observed from the above table that with increase in values of Cb , 

a) N* is increasing. 

b) Mean number of customers in the system is increasing. 

c) Minimum expected cost is increasing. 

Table 10: Effect of C m and C 0 on N*, expected system length and minimum expected cost 



c c 


215,55 


230,60 


245,65 


260,65 


275,70 


290,75 


N* 


9,9 


9,9 


9,9 


9,9 


9,9 


9,9 


L(N*) 


4,4 


4,4 


4,4 


4,4 


4,4 


4,4 


T(N*) 


31.97,3 
4.04 


34.03,34. 
17 


34.08,3 
4.3 


34.13,3 
4.42 


34.19,3 
4.55 


34.24,3 
4.67 



It is observed from the above table that with increase in values of C m and C 0 , 

a) N* is insensitive. 

b) Mean number of customers in the system is insensitive. 

c) Minimum expected cost is increasing. 

Table 11: Effect of Chand C s on N*, expected system length and minimum expected cost 



Ch, Cs 


6,1100 


7,1200 


8,1300 


9,1400 


10,1500 


11,1600 


N* 


8,9 


8,10 


7,10 


7,11 


6,11 


6,11 


L(N*) 


4,4 


3,4 


3,5 


3,5 


3,5 


3,5 


T(N*) 


37.74, 
36.03 


42.2, 
38.06 


44.42, 
40.00 


47.41, 
41.87 


50.20, 
41.68 


52.84, 
45.43 



It is observed from the above table that with the increase in values of Ch and C s 

a) N* is decreasing. 

b) Mean number of customers in the system is decreasing with C h and increasing with Cs. 

c) Minimum expected cost is increasing with both C h and C s . 

VII. CONCLUSIONS 

• Two-phase N-policy M/M/l queueing system with server breakdowns and balking is studied. The closed 
expressions for the steady state distribution of the number of customers in the system when the server is at 
different states are obtained and hence the expected system length is derived. 

• Total expected cost function for the system is formulated and determined the optimal value of the control 
parameter N that minimizes the expected cost. 

• Sensitivity analysis is performed to discuss how the system performance measures can be affected by the 
changes of the both non -monetary and monetary input parameters. 
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